Patterns of Striped Order in the Classical Lattice Coulomb Gas 
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We obtain via Monte Carlo simulations the low temperature charge configurations in the lattice 
Coulomb gas on square lattices for charge filling ratio / in the range 1/3 < / < 1/2. We find a 
simple regularity in the low temperature charge configurations which consist of a suitable periodic 
combination of a few basic striped patterns characterized by the existence of partially filled diagonal 
channels. In general there exist two separate transitions where the lower temperature transition 
(T p ) corresponds to the freezing of charges within the partially filled channels. T p is found to be 
sensitively dependent on / through the charge number density v = pi/qi within the channels. 
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Competition between periodic pinning potential and 
repulsive inter-particle interactions often generate inter- 
esting frustration effects on the equilibrium properties of 
particle systems resulting in diverse commensurate or in- 
commensurate phases In the simplest case of one 
dimensional system on a lattice, there exists a well known 
prescription of obtaining the ground state configuration 
jl],^| for a given particle density provided that the inter- 
particle potential is convex and vanishing at infinity. In 
two dimensions, systematic or analytic results are scarce 

i- 

One of the extensively studied model systems in two 
dimensions is the classical lattice Coulomb gas (LCG), 
where a fraction / = p/q (p and q are relative primes) of 
the lattice sites are occupied by unit charges upon uni- 
form neutralizing background ||. LCG appears in the 
vortex representation of the two dimensional Josephson 
junction arrays under uniform external magnetic field via 
Villain approximation || of the uniformly frustrated XY 
(UFXY) model, where the parameter / corresponds to 
the ratio of the magnetic flux per unit plaquette and the 
superconducting flux quantum <f>o = hc/2e. In these sys- 
tems, physical properties such as transition temperatures 
or critical currents can depend sensitively on (the ratio- 
nality of) the value of the frustration parameter / due to 
the commensuration effect |jj . 

In spite of many years of research effort, there still re- 
mains wide range of / parameter (notably for 1/3 </ < 
1/2) for which not much is known in terms of equilib- 
rium properties. Even the ground state configurations 
are known only for some limited cases of values of /, es- 
pecially for low order rationals [fl3"|| . In this paper, we 
report simulation results on the low temperature ordered 
structures of LCG on square lattices for 1/3 </ < 1/2. 

In the case of UFXY model on square lattices, Halsey 
proposed staircase states (^] as the low temperature or- 
dered configurations with periodicity q x q. However, 



the staircase state turns out to be the true ground state 
configuration only for some limited values of / such as 
/ = 1/2, 1/3, 2/5, 3/8, etc. More recent works on UFXY 
model M and LCG jl0| shows that for / = p/q near 1 — g 
(g = (v5 — 1)/2) with q large, the low temperature vortex 
configuration in the UFXY model can be different from 
the ordered charge configurations in LCG. Especially, in 
the case of LCG on square lattices, it was shown [pLOfl 
that there appear diagonal striped configurations with 
partially filled diagonals. However, the exact patterns of 
low temperature charge configurations for general / of 
dense charge filling have not been enumerated yet. 

In this paper, we find numerically for LCG on square 
lattices with / in a large part of the range 1/3 < / < 1/2, 
that, below a first order transition temperature T c , the 
charge configuration consists of periodic arrangements of 
a few basic striped patterns. Four regimes of values of 
/ are identified, each of which represents a specific class 
of striped patterns characterized by the existence of par- 
tially filled diagonals. T c is found to be only weakly de- 
pendent on /. Below T c , there exist in general a tempera- 
ture range where the charges within partially filled diago- 
nal channels are disordered (thus mobile along channels). 
We found that there exist another transition at a lower 
temperature T p at which freezing of charges within par- 
tially filled channels occur. An interesting result is that 
this lower transition temperature depends sensitively on 
the rationality of the charge number density within chan- 
nels v = pi/qi, decreasing monotonically in q\. 

General 2D LCG [jllj is described by the following 
Hamiltonian, 



i)Qi 



(i) 



where is the distance between the sites i and j, and 
the magnitude of charge Qi at site i can take either 1 — / 
or — /. Charge neutrality condition J^. Qi = implies 
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that the number density of the positive charges is equal 
to /. We can thus view the system as a lattice gas of 
N ■ f charges of unit magnitude upon uniform negative 
background charges of charge density — / (N = L 2 is 
the total size of the system with the linear dimension L). 
The lattice Green's function Gfaj) solves the equation 
(A 2 — \~ 2 )G(rij) = —2n5r zj fl, where A 2 is the discrete 
lattice Laplacian and A is the screening length. For the 
case of usual Villain transformation of UFXY model, we 
have A = oo. But the screening term is included in this 
equation for generality. On a square lattice with periodic 
boundary conditions, G[f) is given by 
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N t—' 2 — cos k x — cos ky + 1/ A 2 ' 



(2) 



where k are the allowed wave vectors with — 2itn ll / L 
(fx = x,y), with = 0, 1, . . . , L — 1. In the case of 
infinite screening length, for large separation r, one gets 
G(r) ~ — lnr JT^j. In this work, the presented results are 
all obtained for the case of A — > oo |l2j . 

In our MC simulations, the initial disordered ran- 
dom configuration is updated according to the standard 
Metropolis algorithm by selecting a positive charge at 
random and moving it over to one of the nearest neigh- 
bor (NN) or next nearest neighbor (NNN) sites |10 . An 



important aspect of our simulations is that one has to 
choose the lattice size appropriately in order to match 
the periodicity of the low temperature configuration (see 
below). If, otherwise, one chooses a lattice size that is 
incommensurate with the periodicity of striped patterns, 
then one ends up with defective charge configurations. 

To begin with, let us present the four basic component 
patterns (Fig. 1) before going into the full description of 
the striped configurations. First component pattern (I) is 
a sequence of three diagonals which are empty, filled, and 
empty respectively (that may be denoted by (010) in our 
notation where 1 refers to a filled diagonal and refers 
to an empty diagonal). In other words, it is a pattern 
with single isolated diagonal filled with charges, that is 
neighbored by empty diagonals on both sides. Repetition 
of this pattern produces the ground state configuration 
for f = 1/3 with spatial periodicity three. 

Second component pattern (II) consists of a sequence 
of five diagonals that can be written as (01010). This 
forms the basis of the ground state configuration for 
/ = 2/5 with lattice periodicity five. The third com- 
ponent pattern (III) consists of a sequence of seven diag- 
onals that can be denoted by (OlOpOlO) where p refers 
to a partially filled diagonal where only part of the di- 
agonal sites are occupied by positive charges. This is 
essentially a partially filled diagonal enveloped by one 
filled diagonal on both sides at second neighbor diagonal 
position, which may be termed as a channel structure. 
This can form a basis of a periodic configuration with 
spatial lattice periodicity seven. Lastly, the fourth com- 



ponent pattern (IV) can be denoted by (OlOlOpOlOlO). 

This is component pattern consisting of a partially filled 
diagonal bounded by two type II stripe patterns on both 
sides. This can form a basis of a periodic configuration 
with spatial lattice periodicity eleven. 

We are now in a position to be able to give a detailed 
description of the low temperature charge patterns. We 
may identify four regimes for the striped charge config- 
uration. We may call these by regime A, B, C, and D 
respectively. Fig. 2 shows typical representative configu- 
ration in each of the four regimes. 

In regime A that is bounded by 1/3 ^ / ^ f c \ 
(fd ~ 0.357), the low temperature charge configuration 
consists of combinations of type I and III patterns, where 
I ( I = 1, 2, 3 ■ • • ) copies of type I patterns in sequence 
followed by a single type III pattern (I'lII) forms a basic 
unit, repetition of which forms the whole charge config- 
uration. We can easily see that the lattice periodicity of 
the ordered stripe configuration is equal to pa — 3Z + 7. 
As I approaches infinity, we recover the case of / =1/3. 
As the value of / increases within regime A, the value 
of I decreases monotonically in step-like manner. There- 
fore, regime A is further divided into sub-regimes each of 
which is characterized by a positive integer I. 

Regime B covers the region f cl <; / <; f c2 (f c2 ~ 
0.381), where the ordered charge configuration simply 
consists of repetitions of type III stripe patterns with the 
resulting lattice periodicity ps = 7. In regime C which is 
bounded by f c % < / < 2/5, the low temperature charge 
configuration consists of combinations of type II and III 
patterns, with m (m = 1, 2, 3 • • •) copies of type II stripe 
patterns in sequence followed by a single type III pattern 
(II m III) forms a basic unit. Here, we can see that the 
lattice periodicity of the ordered stripe configuration is 
equal to pc = 5m + 7. As the value of / increases within 
regime C, the value of m increases monotonically in step- 
like manner. The case of / = 2/5 corresponds to the 
limit of m — ► oo. 

Regime D corresponds to 2/5 < / ^ / C 3 (/ C 3 ~ 
0.425), where the unit period of ordered configuration 
consists of combinations of a type IV pattern plus n 
(n = 0, 1, 2, 3, • • •) repetitions of type II stripe patterns 
that may be denoted by 11™ IV. We see that the lattice 
periodicity of the ordered stripe configuration in regime 
IV is equal to po = 5n + 11. As the value of / increases 
within regime D, the value of n decreases monotonically 
in step-like manner. 

Note that the periodicity in the above refers to the 
periodicity of the filled diagonals only, neglecting the true 
periodicity including the charge configurations within the 
partially filled diagonals. The true spatial periodicity 
of the charge configurations can be many times larger 
than the stripe periodicity since we also have to take into 
account the correlation of charge configurations between 
different partially filled channels. 

For any given value of /, we can easily obtain the 
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filling density v inside the partially filled channel us- 
ing the relations / = (I + 2 + v)/{Zl + 7) (for regime 
A and B), f = (2m + 2 + v)j{hm + 7) (regime C) and 
/ = (2n + 4 + ^)/(5n + ll) (regime D) respectively where 
I, m, n and v are as defined above. As the value of / con- 
tinuously increases within one sub-regime, the system in 
the low temperature stable configuration simply adjusts 
itself by accomodating the extra number of charges into 
the partially filled diagonal channels and thereby chang- 
ing the charge filling v within the channels. 

We find that there exists in general another transition 
at lower temperature T p corresponding to charge freezing 
inside partially filled channels. In order to check that, 
we calculated the inverse dielectric constant along two 
perpendicular diagonals to identify the transition tem- 
peratures. The wave-vector dependent inverse dielectric 
constant is defined as follows G3], 

e- 1 W=(l-^<^- fc >), (3) 

where pf. = J\ Q(rj) exp(— ik ■ fl) is the Fourier compo- 
nent at wave-vector k of the charge density and fi = L 2 
is the total area of the system. By letting k — > for 
a given direction of the wave-vector, one can obtain the 
long wavelength limit of the inverse dielectric constant 
along a specific direction. 

Figure 3 shows the dependence of the inverse dielectric 
constant along parallel and perpendicular to the stripes 
respectively for / = 13/35. We can clearly see that there 
exists, in addition to the transition at T c ~ 0.32 cor- 
responding to the onset of striped order, intermediate 
regime of temperature where the dielectric constant ex- 
hibits asymmetry due to the channel-striped structures. 
Also, the system is seen to undergo another transition 
at lower temperature T p ~ 0.014 (determined rather ar- 
bitrarily as the temperature where the inverse dielectric 
constant is equal to 0.4). Figure 4 shows the dependence 
of the inverse dielectric constant along parallel to the 
stripes for various values of / in regime B, where a wide 
variation is seen in the temperature dependence of the 
inverse dielectric constant parallel to the channels. 

Shown in Fig. 5 is the dependence of the two transi- 
tion temperatures on / = p/q in regime B where the 
higher transition T c is seen to depend on the values of 
/ smoothly, while the lower transition temperature T p 
exhibits sensitive dependence on /. From the inset of 
Fig. 5 where a plot of T p versus the integer denomina- 
tor qx of the charge number density v = p\/q\ within 
channels, one can recognize that T p decreases monoton- 
ically as q\ increases. This is another commensuration 
effect coming from the rationality of the particle number 
density within channels each of which forms effectively a 
one-dimensional lattice gas system if we neglect the in- 
teraction between different channels. The arrangement 
of charges within each channel was found to follow the 



pattern given in ref. JjJ and Q at least for the values of 
/ considered in this work. 

One can ask what determines the boundary value of / 
and v for each of the sub-regimes, in other words, what is 
the stability criterion for each pattern of the striped con- 
figuration. Even though an analytic formula cannot be 
given, our simulations suggest that the criterion of de- 
termining the crossover point between two neighboring 
subregimes is related to the electrostatic stability which 
is determined by the filling density inside the partially 
filled channel. We found numerically that for any given 
value of / the filling density v in the channel approx- 
imately satisfies the inequality 0.4 Ss v ^ 0.7. For a 
given sub-regime, possible values of v were always within 
this bound. Beyond some threshold value v c of v that is 
within the above bound, electrostatic instability begins 
to set in, and rearrangement of the whole charge con- 
figuration occurs in order to form a new stable ordered 
patterns. 

Now we are left with the region of / C 3 < / < 1/2. 
Even though we have not investigated the ordered con- 
figurations extensively for all values of / in this regime, 
we could see, from some of our annealing simulations for 
rational values of / in this regime, that the low tem- 
perature configuration no longer shows striped patterns 
but rather consists of regular arrays of hole defects upon 
/ = 1/2 checkerboard configurations. If we suppose that 
/ = 1/2 — a = (1 — pijqi)j2 with p 2 and q 2 relative 
primes, then we can easily see that /' = p 2 /q 2 = 2a 
represents the density of hole defects upon the / = 1/2 
checkerboard configuration. Therefore, we can suspect 
that the low temperature defect configuration in this 
regime will be equivalent to the charge configuration of 
lattice Coulomb gas with / = /'. We could confirm this 
expectation for a few cases of / [ p~3| , p~5|| - 

In summary, we have shown numerically that the 2D 
LCG on a square lattice exhibits a simple regularity in 
its striped charge configuration at low temperatures for 
filling factor / in a large part of the range l/3</<l/2 
which is characterized by the existence of partially filled 
diagonals. The low temperature ordered configuration 
consists of a simple combinations of four basic striped 
patterns. In general, there exists another transition at 
a lower temperature T p corresponding to the freezing 
within partially filled channels. It would be interesting 
to observe these striped charge patterns experimentally, 
e.g., in regular square arrays of ultrasmall tunnel junc- 
tions. It may also be possible to observe similar patterns 
in the macroscopic systems of dielectric charged spheres 
]l(| under periodic pinning potentials. 
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FIGURE CAPTIONS 

Fig. 1. Regimes of charge patterns for the range of value of 
/ between 1/3 and 0.425. Filled squares and empty 
squares represent positive (Q = 1 — /) and negative 
(Q = — f) charges respectively, while gray squares 
denote lattice sites forming partially filled diagonal 
channels, where only finite fraction v of the sites 
are filled with positive charges. 

Fig. 2. Low temperature charge configurations for (a) 
regime A with / = 7/20 (I = 1), (b) regime B with 
/ = 13/35, (c) regime C with / = 0.384 (to = 1) 
and (d) regime D with / = 0.41 (n = 1), respec- 
tively. 

Fig. 3. Inverse dielectric constants along parallel and per- 
pendicular to the stripes for / = 13/35 and L = 
35 versus temperature. We can clearly see an 
anisotropy of the inverse dielectric constant. 



Fig. 4. Inverse dielectric constant along parallel to the 
stripes for various vaues of charge number density 
v within channels. All of the systems shown arc 
chosen from regime B. 

Fig. 5. Transition temperatures (both T c and T p ) versus 
the charge filling ratio /. Inset shows the depen- 
dence of T p on the denominator q\ of v = pi/qi . 



4 




Lee, et al 




Lee, et al 




Lee, et al 



(d) 



Lee, et al 



Lee, et al 




Lee, et al 




0.06 



0.05 



0.04 



0.03 



0.02 



0.01 { M T p vs. f 

A A T vs. f 





0.35 0.36 



0.016 



OO T vs. q, 




0.37 0.38 0.39 

f 



